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1 Introduction 



At zero temperature, the spectrum of QCD is encoded in hadronic spectral func- 
tions. Groundstates, excited states, decay widths and continuum contributions 
can, in principle, be extracted from spectral functions in different channels. Sim- 
ilarly, at finite temperature and/or density, medium modification of hadrons, the 
rate of photon and dilepton production, and hydrodynamical response functions 
can be obtained from the appropriate spectral functions. 

Spectral functions are inherently real-time correlation functions and therefore 
difficult to obtain using standard lattice QCD data analysis techniques. Since eu- 
clidean lattice correlation functions are determined numerically on a finite num- 
ber of points in imaginary time only, the analytical continuation to real time is 
classified as an ill-posed problem. In the past few years significant progress in 
the extraction of spectral functions from lattice QCD has come from the applica- 
tion of the Maximum Entropy Method (MEM) to this problem. MEM has been 
applied in many branches of science (see e.g. Ref. [1] for a review), a thorough 
review focussing on lattice QCD can be found in Ref. [2]. 

In order to interpret hadronic spectral functions obtained in lattice QCD, it is 
important to understand how lattice artefacts will appear. Lattice artefacts are 
expected to be present at large frequencies and may be studied in perturbation 
theory. Free lattice meson spectral functions have been studied at zero momen- 
tum for Wilson and hypercube fermions [3] and at nonzero momentum for Wilson 
and staggered fermions [3]. In this paper our first goal is to study lattice me- 
son spectral functions for chirally symmetric fermions, specifically overlap [HIE], 
domain wall [3, El [9], and overlap hypercube [TQj, [TT] fermions. Furthermore we 
present results for meson spectral functions obtained with the Maximum Entropy 
Method in QCD with 2 + 1 flavours of domain wall fermions, using data generated 
by the UKQCD and RBC collaborations on QCDOC machines. 

The paper is organized as follows. In the next section we derive general ex- 
pressions for free meson spectral functions on a finite lattice, independent of the 
particular fermion formulation that is used. In Section [3] we specialize to overlap 
fermions and compare the resulting spectral functions, and in particular lattice 
artefacts, with continuum and staggered spectral functions. This analysis is ex- 
tended to domain wall fermions in Section HJ and to overlap hypercube fermions in 
Section O In Section E] we present first results for pseudoscalar spectral functions 
in QCD with dynamical domain wall fermions. We find good agreement between 
the MEM results and the groundstate mass obtained with conventional cosh fits. 
We also argue that the structure seen at larger energies is consistent with lat- 
tice artefacts found in the first part of the paper, but a quantitative comparison 
requires further study. 
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2 Spectral functions 



We start with a brief summary of well-known relations |2J. Euclidean meson 
correlators are defined by 

G? ff (r > x) = (J H (r,x)4(0,0)), (2.1) 

where Jh(t, x) = ^(r, x)IV0(t, x) and = {1, 75, 7", 7^75} for the scalar, 
pseudoscalar, vector and axial vector channels respectively. In general euclidean 
(and other) correlation functions are related to their spectral function via a dis- 
persion relation in momentum space, 

Gh{z ' p) -J 2^^7^' (2 - 2) 

where z is the frequency extended into the complex plane. Equating z to iu n , 
where u n = 2irnT (n G Z) is the Matsubara frequency, yields the euclidean 
correlator at finite temperature T. In imaginary time this relation reads 

G h (t, p) = J™ ^ K(t, u)p H (u, p), (2.3) 

with the kernel 

cosher -1/2T)] 

K{t,u) = . , „>. • (2.4) 

sinh(u;/2T) 

At zero temperature, this kernel reduces to K(t, uj) = e~ WT . 

In this section we derive a general expression for meson spectral functions 
on an isotropic lattice with N% x N T sites, by writing the euclidean correla- 
tors in the form (I2.3P and identifying the lattice spectral function from that 
expression. We use periodic boundary conditions in space, fcj = 2ixni/N a with 
Hi = -NJ2 + 1, -N a /2 + 2, . . . , NJ2 - 1, N a /2 for 2 = 1,2, 3, and antiperi- 
odic boundary conditions in imaginary time, k^ = 7r(2n4 + 1)/N T with n.4 = 
-N T /2 + 1, -N T /2 + 2,..., N T /2 - 1, N T /2. Lattice units a = 1 are used through- 
out. 

The correlators we are interested in are of the form 

Gh(t, p) = ~ tr S ^ k ) T HS(-r, p + k)T H , (2.5) 
a k 

where S(t, k) is the fermion propagator and N c denotes the number of colours. In 
order to derive compact expressions for lattice meson spectral functions, starting 
from Eq. (12. 5p . we first discuss a generic fermion propagator, without specifying 
a particular formulation. We consider the following fermion propagator 



S(K) 



D{K) 



4 

— 1 

v=l 



J2 c A K )lu sin k v + m{K) 



(2.6) 
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where K denotes the four-momentum. The functions C V (K), m(K) and D(K) 
depend on the fermion formulation, but they are all invariant under k u — > — k v . 
In order to be able to use Eq. (12. 51) . we construct the fermion propagator in the 
mixed representation, 

S(r,k) = ±-J2 elk4TS ( K )- ( 2 - 7 ) 

We assume that S(K) has a single pole at fc 4 = HE^, determined by D(iEy_, k) = 
0. In the case of more than one pole, a summation over the poles is required. 
This yields [H] 

3 

S(t, k) = 74 S 4 (r, k) + iMt, k) + 1S u (t, k), (2.8) 

i=l 

where 

Si{r, k) = >S 4 (k) cosh(f£'k)? 
Si(r, k) = Si(k) sinh(f£' k ), 

S u (r,k) = S u (k)sinh(f£ k ). (2.9) 

Here < t < N T = 1/T and f = r — 1/2T. The momentum-dependent coeffi- 
cients read 

C 4 (iE k , k) sinh£ k 

4 (KJ — 



Si(k) 



2£ k cosh(£ k /2T) ' 
C^iE^, k) isinfcj 
2^ cosh(£ k /2T)' 



where 



S„(k) = - "^Ly (2.10) 
2£ k cosh(i^ k /2i ) 



11 , 

i Res — — -. (2.11) 



2£ k ki=iE k D(K)' 

The propagator satisfies 5(— r, k) = 7sS''''(r, k)7 5 . 

For reference, we note that for fermions in the continuum one finds the 
same propagator, with the replacements £ k — > _E k , sinhi*^ — > i? k , sin 
m(iE k , k) — ► m and — > 1. 

Inserting Eq. (12.81) in Eq. (12.51) gives the euclidean correlator 



GH(r,p) = -^^[a«^ 4 (r,k)5i(r,r) 
CT k 

-ag } S<(t, k)5j(r, r) - ag^r, k)Sj(r, r)l , (2.12) 
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a (1) 


(2) 
a H 
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a (1) 


(2) 


(3) 
% 


PS 


1 


1 


-1 


1 


Pps 


75 


1 


-1 


-1 


p oo 


7° 


1 


1 


1 


P^° 


7°7s 


1 


1 


-1 


p ii 


y 


3 


-1 


-3 


Pf 


7*7s 


3 


-1 


3 


Pv 




2 


-2 


-4 


Pa 


7^75 


2 


-2 
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Table 1: Coefficients a H for free spectral functions in different channels H. In 
the case of Y and 7*75, the sum is taken over i = 1,2,3; pv = p w — P°° and 

„ „ii „00 

PA ~ P 5 - P 5 • 



where r = p + k. The coefficients ag* are given in Table [TlFl 

We will now extract the lattice meson spectral functions. It is straightforward 
to write the above expression for Gh(t, p) as 



Gh(t, P )= r d ^K{r^)p 

Jo 27T 



lattice/ \ 



(2.13) 



and identify the expressions for the lattice spectral functions [I] 



rif'"(-,P) = ^£ S inh V2r , 

a k 



+ 



g^k^r) + ag> £ S,(k)5}(r) + ag^k)^ (r 
ag^k^r) - ag } ^5i(k)5j(r) - ag^k^r) 



<J(u; + .Bk - E T ) 
8{u -E^-E r 



-OJ. 



(2.14) 



This expression is the immediate counterpart of the continuum result [4] 

d 3 k 1 



p c ™ t (cu,p) = 27rN c J 



(2tt) 3 E k E r 

[n F (E k ) - n F (E r )} ag *E k E r + ag } k ■ r + agW 5{uj + £ k - E r ) 

+ [l-n F {E k ) -n F (E r )} a% ] ' E k E r - ag } k • r - ag } m 2 5{u - E k - E r ) 

-u)\, (2.15) 



(CO 



where n F (uj) = l/[exp(u/T) + 1] is the Fermi distribution, as can be seen by 
making the appropriate substitutions. At zero temperature, only the "1" in the 



1 We use Minkowski gamma-matrices to label the channels, see Ref. [4] for further details. 
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second term survives. In the continuum the remaining three-dimensional integral 
can be carried out for arbitrary external momentum and quark mass. Analytical 
expressions for continuum meson spectral functions can be found in Ref. [4]. 
Since the expressions are lengthy, we give here the result at vanishing external 
momentum only, 



P r(^0) = e(^-4m 2 )^/l- 



in 



lu a 



< u ay 

-Anuj5(uj)N c 



im 2 ( a% ] 



1 ■-■2n,,[ - 



d 3 k 



n F (u k ) 



tot. 



+ m 



,(!) 



H + a H 

(2.16) 



The first term contributes above threshold {u > 2m), while the contribution pro- 
portional to uj8(uj) is related to conserved quantities (see e.g. Ref. [15] in relation 
to transport coefficients). Note that spectral functions are odd, p H (—u,p) = 
-p H (uj,p). 

On the lattice the spectral functions (12.141) can in general not be evaluated 
analytically. Instead, in the following sections we give the explicit expressions for 
the free fermion dispersion relation E^, the coefficients S^k), Si(k) and 5 M (k), 
and the residue S^, for overlap, domain wall and overlap hypercube fermions. We 
use those in Eq. (12.141) to compute lattice spectral functions by performing the 
spatial lattice sum over k numerically, using the same approach as in Refs. [31 H]. 

3 Overlap fermions 

The massless overlap (Neuberger) operator is given by 



D N = n [ 1 + ^==) ' X = D w - fi, 



(3.1) 



where Dw the usual Wilson-Dirac operator and [i is a constant]^ In momentum 
space, X reads 

X{K) =i^ lv smk v + b{K), b(K) =rj^(l- cos k v ) - M- (3-2) 



The corresponding fermion propagator is 



s(K) = -1 ( - l h lu Z k : + 1 



2/i V. oj(K) + b(K) 

2 In other studies the coefficient /i is sometimes denoted by p 



(3.3) 
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with 

u{K) = VX^X = l^2sm 2 K + b 2 (K). (3.4) 

The poles of the propagator are determined by ou(K) + b(K) = 0, which yields 

^sin 2 ^ = 0, b(K)<0. (3.5) 

V 

Writing = iE^ gives the dispersion relation 



coshE k = yjl + /C 2 , b{iE k ,k)<0, (3.6) 

where we defined 

Kl = J^sin 2 ^. (3.7) 

i 

The constraint b < arises from the square root in the definition of Due 
to this constraint, there is only a pole when k is not too large, and there are no 
solutions near the edges of the Brillouin zone. Apart from this, the dispersion 
relation is identical to the one for naive fermions. 

Meson spectral functions for free massless overlap fermions take the form 
(I2.14p of the previous section, with the now explicitly determined functions (12.101) 

a(^k,k) = l, S u {k)=0, (3.8) 

and the residue (12.111) 

£k /i sinh E k cosh E k 

Chiral symmetry is manifest in the meson spectral functions, since the expressions 
in the scalar (vector) and the pseudoscalar (axial vector) channel only differ with 
respect to the coefficient Table [TJ Since S u = 0, dependence on 

vanishes for massless overlap fermions. 

The extension to massive overlap fermions is straightforward. To include the 
mass, the overlap operator is changed to 

( m o\ 

D O v,m = 1 - — ) D N + m , (3.10) 



2/i 

and the corresponding propagator is 



c ( K) _ 1 jfj ~ m °/ 2 ) H Eg 7. sin K + b(K)} + {n + m /2)u(K) 

2 {v? + ml/A)u{K) + {^-ml/A)b{K) ' ' 



3 Provided that < /i < 2r, there is no pole at = ir — iE^, since then the constraint 
cannot be met. 
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We find that the pole is determined by 



J^sin 2 k v = -m 2 b 2 (K), 

V 

where we defined 



{fi 2 - ml/ '4) b(K) < 0, (3.12) 



m 



fim 



/i 2 + mg/4- (3 - 13) 

We always consider tuq to be small (in lattice units), while \x ~ 1, so that the 
constraint still implies b < 0. Writing again = iE k and solving the resulting 
quadratic equation gives the allowed dispersion relation 



cosh Ek 



1 — r 2 mg |_ 



rrh 2 Mk + J (1 + /C 2 )(l - r 2 m 2 Q ) + rh 2 Ml 



provided that 



Here we have defined 



b(iE k , k) = -r coshE k + M k < 0. 



A^k = r — fi + r — cos ^0- 

i 

At zero momentum, the rest mass is determined by 

1 



(3.14) 
(3.15) 
(3.16) 



cosh Eq 



1 — r 2 m,Q 



r(fi — r)rriQ + yl + mlnin — 2r) 



while at small k and m Eq. (I3.14p reduces to 

coshEk = 1 + - (k 2 + m 2 ) + . . . 



(3.17) 



(3.18) 



as expected. 

We compare the massive overlap dispersion relation with the naive (or stag- 
gered) dispersion relation, 



cosh E k = J 1 + K. 2 , + ml, 



(3.19) 



and the continuum expression E k = ^k 2 + ml in Fig. [T] for two values of the 
overlap parameter /i and a rest mass ttir = Eq = 0.1 (the HF kernel shown as 
well will be discussed below). Throughout this paper we take r = 1. We find 
that the overlap dispersion relation can hardly be distinguished from the naive 
one for the value of ttir shown here. In the overlap case the dispersion relation 
terminates before the edge of the Brillouin zone due to the constraint b < 0. 
The momentum value of the endpoint depends on the overlap parameter fi and 
increases with increasing fi. 
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Figure 1: Dispersion relation along three directions in the Brillouin zone for 
massive overlap fermions with rest mass ttir = 0.1 using a massive/massless HF 
kernel with and a standard Wilson kernel with fi = 1, 1.8. For comparison 

the naive and continuum dispersion relation are shown as well. 
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0.5 1 1.5 2 2.5 



Figure 2: Pseudoscalar spectral functions pps(w, 0)/u; 2 for staggered and standard 
overlap fermions with \i = 1, 1.8, and m R = 0.1, N T = 32. 



The coefficients in the meson spectral functions now read 

C„(i£ k , k) = 1 - ^, m(iE k , k) = m ° ,M iE*, k), (3.20) 

and the residue is 

1 n 1 u(iE k ,k) 



£\l /i 2 + mp/4 cosh E k + rriQrb(iE k , k) sinh 



(3.21) 



Comparison between these functions at small momentum and their continuum 
counterparts shows that there is a multiplicative renormalization of the fermion 
propagator at finite itlq. This is not unexpected since the kinetic term in (13.111) 
has a nonstandard normalization. One way to write the renormalization factor is 
to compare the expressions at zero spatial momentum. Explicitly, the renormal- 
ization factor is given by 

C v (iE ,0)^-. (3.22) 

At small K and mo, the euclidean fermion propagator (13.111) is approximately 
given by 



S{K) « ^~ mo/2 2 } 
(/i 2 + m 2 /4) 2 



-iJ/C + m 1 



K 2 + m 2 2/x 



(3.23) 



The multiplicative prefactor goes to 1 in both the chiral and the continuum limit. 
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a(0 



Figure 3: N T dependence of the pseudoscalar spectral function pps(u;, 0)/u 2 for 
standard overlap fermions (/x = 1, m R = 0.1). 



We now have all the ingredients to compute spectral functions with free over- 
lap fermions. In Fig. [2] we show the pseudoscalar spectral function at zero mo- 
mentum for overlap fermions with /i = 1 and 1.8, staggered fermions and the 
continuum result (12.161) . All spectral functions increase in the same manner be- 
yond the threshold u = 2m^, provided that tur is small. At larger frequencies, 
effects due to the deviation of the continuum and lattice dispersion relation be- 
come visible. The first cusp is due to the maximal lattice energy reached at 
k = (7r/2,0,0) (cf. Fig. [D top), which yields a cusp at au « 2 cosh" 1 y/2 « 1.76. 
The difference in height of the spectral functions is due to the different residues, 
which depends on \x for overlap fermions and for staggered fermions is given by 

m 

£k cosh Ek sinh ^ 

For staggered fermions, there is a second cusp due to the maximal lattice energy 
reached at k = (7r/2, tt/2, 0). For overlap fermions this cusp is absent (for \i = 
1) or less pronounced (for \x = 1.8), since the constraint has terminated the 
dispersion relation (cf. Fig. [1] middle). The maximal energy is given by co> max = 
2-Ekmax, which for staggered fermions is reached at k = (n/2, tt/2, tt/2) and given 
by auj w 2 cosh" 1 2 m 2.63. For overlap fermions, this region is again modified 
due to the constraint. 

The effect of finite iV T is shown in Fig. [3] for overlap fermions with \i — 1 and 
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is seen to be negligible for large enough values of N T . 

Finally, lattice spectral functions at nonzero momentum are similar to the 
continuum ones in the small frequency region au < 0.5 and lattice artefacts at 
larger w are not affected by the external momentum, provided that the external 
momentum is small (see Ref. 0] for further details). 

4 Domain wall fermions 

The four-dimensional massive domain wall propagator, assuming an infinite ex- 
tent of the fifth dimension, reads 





u sin k v + m (1 




W\ 


e~ a ) 


e a \ 


W\ 


-l + m§(l- 


W 


\e~ 


■a) 



where 

l + W 2 + Y^ sin 2 k ( . 
cosh a = — , (4.2) 

W = l-// + ^(l-cosfc„)- (4.3) 

V 

The domain wall height is denoted with \x and is taken between and 2. For 
small momentum and 1 < fi < 2, W is negative (see e.g. Ref. [16] for a recent 
review) . 

The dispersion relation is determined by the pole in the propagator at k^ = 
iE k . Again we find a quadratic equation for coshi*^, with the allowed solution 

coshE k = Xk + (1 + m ^ , (4.4) 

Z k 

with 

x k = -2m 2 (1 + M k ) [1 + Kl + M k {2 + M k )] , 

y k = (1 + /C 2 ) (1 - ml) 2 + ml [l + /C 2 - M k {2 + M k )] 2 , 

z k = (l- m 2 ) 2 - Am 2 M k (2 + M k ), (4.5) 

where M. k is given in Eq. (13.161) with r = 1. 

For massless domain wall fermions, Eq. (I4.4p reduces to 



cosh£ k = ^/l + ZC 2 , (4.6) 
as in the overlap formalism. The rest mass is determined by 



2m 2 (/i - 2) 3 + (1 + m 2 ) + m 2 (2 + gjj - 4)Qu - 2)") + m\ 
C ° ShE ° = l-2m 2 (2^- 8 ^ + 7)+m^ " { ^ ] 
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Expanding the massive case for small k and mo, we find 

cosh£ k = l + ^(k 2 + m 2 ff ) +..., (4.8) 

where m e fr = (1 — wfymo is the multiplicatively renormalized fermion mass. The 
multiplicative factor reads 

l-wl =/i(2-/i), w = W{0) = l-fi. (4.9) 

In the limit that tuq <C 1, tur = m e ff. The constraint on the allowed fermion 
energies arises in this case from the behaviour of the propagator in the fifth direc- 
tion: demanding a normalizable solution of the five-dimensional Dirac equation 
yields the constraint 

\W(iE k ,k)\<l, ( 4 -!0) 

both in the massless and the massive case0 

The resulting dispersion relation ( 14.41) is indistinguishable from the dispersion 
relation for massive overlap fermions for small values of mo, shown in Fig. [TJ At 
finite values of mo, deviations are more pronounced for larger values of \i. 

The coefficients in the meson spectral functions read 

C v (iE k , k) = 1, m(iE k , k) = mo (l - , (4.11) 

and the residue reads 

1 _ 2 VA 2 - AW 2 

£k ~ A + (1 + M k ) -2W + m 2 \Al (1 + M k ) - 2W] sinh E k 

Here we defined 



(4.12) 



A± = A±VA 2 -4W 2 , A = 1 + W 2 - sinh 2 £ k + /C k , (4.13) 

and all quantities are evaluated onshell at K = (iE^, k). In the massless case the 
residue simplifies considerably to 

1 = i-w 2 m,k) ^ 

£ k cosh E k sinh E k 
At small K and m , the euclidean fermion propagator is approximately given by 



4 We note here that W(i£'k,k) is always larger than — 1, so that this constraint coincides 
with 6(i-Ek, k) < in the overlap formalism, since W = 1 + b. 
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a(0 

Figure 4: Pseudoscalar spectral functions pps(w, 0)/u; 2 for domain wall fermions 
(with /i = 1, 1.3, 1.4 and // = 1.8 in the inset) and standard overlap fermions 
(with jj, — 1). 



which suggests a multiplicative renormalization factor 1 — Wq. However, we have 
found numerically that corrections to this renormalization factor are large for 
small but finite m and larger values of \l. From a comparison of the expressions 
at zero spatial momentum, we find the multiplicative factor to be 

f 2 , (4.16) 

which deviates substantially from 1 — w\ for the largest value of \i used here. 
The resulting spectral functions are shown in Fig. |4] for different values of /i. To 
obtain these spectral functions, we fix the rest mass m# = 0.1 and solve for the 
bare mass m using Eq. ( 14. 7ft . Subsequently we find the multiplicative factor 
using (I4.16P to properly set the vertical scale. 

As can be seen in Figs. H] and [2j the spectral functions obtained from overlap 
and domain wall fermions differ at larger frequencies. The overlap spectral func- 
tions show reduced discretization effects. We note that is due to the difference in 
residues, and not because of the dispersion relations. Increasing the domain wall 
height shows that the \x dependence is much stronger than in the overlap case. 
For /i = 1.8 the effect is remarkably large. 
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5 Overlap hypercube fermions 



The overlap formalism solves the chirality problem in lattice QCD. However, its 
dispersion relation shows no reduction in discretisation effects when compared to 
e.g. staggered fermions. A systematic approach to constructing fermion actions 
which have good chiral properties and very small discretisation errors is based 
on the renormalisation group (RG). In this approach one aims to approximate 
so-called 'perfect actions', which lie on the renormalized trajectory of an RG 
transformation [T71 [18]. Here we consider a particular truncation of a perfect 
action, the hypercube fermion (HF) action introduced in Ref. [19]. This action 
is obtained from an RG transformation which yields an ultralocal action in one 
dimension [20] . For a recent review, see Ref. [21 J . 
The hypercube action is written as 

s = ^ip(x) I ^ivpv(x -y) + K x - y) \ *l>{y)> (5- 1 ) 

x,y {. v ) 

where the couplings p v and A are nonzero only if x and y are in the same hyper- 
cube. Explicitly, the corresponding Dirac operator reads, in momentum space, 

D RF (K) = tJ2C? F (K) lu smk u + m HF (K) } (5.2) 

with 

C* F (K) = 2p (1) +4p (2) ^cosA; (T + 8p (3) ^ JJ cos k v + 16p (4) JJ cos k a , (5.3) 
and 

m RF (K) = A (0) + 2A (1) cos k " + 4A(2) II cos k *> 

+8A (3) n c ° s k ° + i6A(4) n c ° s k »- ( 5 - 4 ) 

The coefficients and A^ a ^ are determined by requiring that this action repro- 
duces the perfect action on a finite volume with sides of length 3, and periodic 
boundary conditions, see Table [2] for two examples. 

In the limit of zero momentum, the action coefficients satisfy 

CHO) = 2 [p« + 6p( 2 ) + 12P® + 8p( 4 >] = (^t) 2 , (5.5) 
mHF (0) = A<°> + 8A« + 24A® + 32A< 3 > + 16A< 4 > = (5.6) 
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massless HF 


massive HF 




massless HF 


massive HF 


p(l) 


0.136846794 


0.054580 


A (o) 


1.852720547 


1.268851 


p (2) 


0.032077284 


0.011010 


AW 


-0.060757866 


-0.030083 


p (3) 


0.011058131 


0.003255 


A (2) 


-0.030036032 


-0.010830 


p (4) 


0.004748991 


0.001206 


A (3) 


-0.015967620 


-0.004716 








A (4) 


-0.008426812 


-0.002212 



Table 2: Coefficients and in the hypercube action for massless = 0) 
and massive = 1) hypercube fermions. 



where tjir = Eq is the rest mass in HF dispersion relation. These relations follow 
from the expression for the one-dimensional fixed point action, which can be 
evaluated explicitly, and ultimately they depend on the RG transformation used 
to construct the action. 

The truncation involved in the construction of the hypercube action intro- 
duces chiral symmetry breaking and discretisation errors. If the truncation is 
justified, these effects will be smallj^l Following Ref. [11] , exact chiral symmetry 
can be restored by using the hypercube operator as the kernel for the overlap 
operator. The resulting overlap operator should inherit many of the properties of 
the kernel and, in particular, have much smaller cutoff effects than the standard 
overlap operator. 

To determine the expression for the overlap hypercube propagator we first 
write the expression for the kernel, X = Dhf — in momentum space, and 

X(K) = i J2C* F (K) liy sin k u + b(K), b(K) = m UF (K) - //. (5.7) 

The corresponding propagator for massive overlap fermions is obtained by mul- 
tiplying sinfc„ with C^ F (K) in Eq. (IBTTTj) of Section H Writing k A = iE kl the 
dispersion relation is again determined by a quadratic equation for cosh E^. Since 
the explicit expressions are rather lengthy, we have listed them in Appendix [A] 
In the limit of small k and mo, we find that the overlap HF dispersion relation 
reduces to 

cosh E k = 1 + - (k 2 + m 2 ff ) + . . . (5.8) 

where the renormalized fermion mass is 

/i-m HF (0) 

meff = KT(o) m °" ( } 

Using Eqs. (15.51 15. 6p . we find that the overlap mass m receives a multiplicative 
renormalization for massive HF, while for massless HF such renormalization is 



5 For studies of meson spectral functions using hypercube fermions in quenched QCD, see 
Refs. [TJ [13]. 
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Figure 5: Pseudoscalar spectral functions pps(^, 0)/u; 2 for overlap fermions with 
a standard and a massless/massive HF kernel (fi = 1, itlr = 0.1, N T = 32). 
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Figure 6: As in Fig. [5] for the vector spectral function py(ui, 0)/u 2 . 
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absent. For massless overlap fermions, the dependence on the coefficients in the 
HF kernel cancels completely in the limit of small momentum. 

The resulting dispersion relations along three direction in the Brillouin zone 
are shown in Fig. [1] for overlap fermions with both a massless and a massive (with 
HF rest mass rnn = 1) HF kernel. Note that we fixed the overlap rest mass at 0.1 
and determined the bare mass from Eq. (I5.9P . It is clear that the improved scaling 
of the HF kernel ensures agreement with the continuum dispersion relation for 
much larger momenta. The deviation at larger momenta results in the behaviour 
dEk/dk > 1, which is especially pronounced for the massless kernel. 

The coefficients in the meson spectral functions are given in Eq. (I3.20p . af- 
ter multiplying C^iE^, k) with C^ F (iE k , k). The residue is given by Eq. (1A.9j) . 
Comparison of the coefficients and the residue with their continuum counterparts 
shows again that the fermion propagator receives a multiplicative renormaliza- 
tion. Expanding the HF overlap propagator for small K and m yields 



S{K) 



/i 3 (/i - Wp/2) 



+ 



(5.10) 



where m e s was defined above. For overlap fermions with a massive HF kernel 
we find therefore that both the overlap mass and the fermion propagator are 
renormalized by the factor \p — wt-hf (0)] //xC^ F (0). 

Meson spectral functions with overlap hypercube fermions are shown in Fig. [5] 
for the pseudoscalar and Fig. [6] for the vector channel. Besides a factor of 2 (see 
Tabled]), these spectral functions differ in detail around the knee at auo ~ 0.5, 
which can be understood from the continuum expression (12.161) . As expected, the 
first cusp in the spectral functions is shifted to substantially larger frequencies, 
determined by twice the maximal energy along the (1,0,0) direction, yielding 
auj ~ 3.75 (5.0) for the massive (massless) HF kernel, as can be seen from Fig. [1] 
(top). As a result, the continuum behaviour at larger frequencies, p{oj) ~ u 2 , is 
better reproduced. However, we would like to point out that the contributions 
from these large frequencies are highly suppressed in the euclidean correlator. 
Taking for simplicity the zero temperature kernel K(t,u) = e~ UT = e~ auJ1lT , we 
find that already at the first time slice K ~ e -5 when au ~ 5, demonstrating the 
insensitivity to these large frequencies. 



6 QCD with dynamical domain wall fermions 

We now build on the free field calculations and consider spectral functions in 
QCD. Most spectral function studies to date have been carried out in quenched 
QCD; for recent work, see e.g. Ref. [22J and references therein. A study of charmo- 
nium spectral functions in dynamical QCD with two flavours on highly anisotropic 
lattices can be found in Ref. [23] ■ Spectral functions at nonzero momentum in 
quenched QCD are considered in Refs. [2"H [25]. 
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In this section, we consider spectral functions at zero momentum, extracted 
from meson correlators obtained in dynamical lattice simulations. These corre- 
lators are computed using the domain wall fermion action on 2+1 flavour back- 
ground configurations employing the Iwasaki gauge action [26], [27] J§ This data 
has been generated by the RBC and UKQCD collaborations [2SI EH EH E21 E3] 
on QCDOC [311 [35], [36]. Here, we present results obtained on a 16 3 x 32 lattice 
at (3 = 2.13 with an inverse lattice spacing of 1.6 GeV [33J. The number of points 
in the fifth dimension is N s = 16. 

In practical implementations of the domain wall fermion formalism it is im- 
portant to choose a domain wall height that minimises the mixing between right 
and left-handed fermion modes at zero bare quark mass. This mixing, which 
vanishes in the limit iV s — > oo, induces some residual chiral symmetry breaking 
and generates, for example, an additive quark mass renormalisation. In free field 
theory the optimal value for the domain wall height is unity. However, in the 
interacting theory, the domain wall height receives a large additive correction and 
the bare parameter should be adjusted accordingly (37J [38]. The bare domain 
wall height used in these simulations is /i = 1.8. Subtracting a simple mean-field 
estimate for the radiative corrections [39] yields a value for the domain wall height 
of /x MF = 1.3029 [40], which is much closer to unity. 

In this first study we show results obtained with a light bare sea quark mass 
wind = 0.02 and a heavier sea quark mass m s = 0.04. To preserve unitarity, 
the valence quarks are constrained to take the same bare mass values as the sea 
quarks. To determine the spectral functions the Maximum Entropy Method [2] is 
applied to correlation functions measured on an ensemble of 70 independent gauge 
field configurations. The meson interpolating operators used are local quark field 
bilinears. The meson correlators are symmetric about the central time slice of 
the lattice. In our analysis we average the correlation functions over equivalent 
time slices and exclude the contact term at r = 0. The MEM algorithm uses 
Bryan's method [H]. As the default model we take Pdefauit(^, 0) = m uj 2 , where 
the constant m is determined by a best fit to the data. 

In Fig. [7J we show spectral functions obtained from pseudoscalar correlators 
evaluated for different valence quark masses. For each quark mass combination 
clear peaks are visible whose position corresponds to the energy of the lightest 
state that couples to the interpolating operators. The horizontal errorbars on 
these peaks give a measure of the peak width, while the vertical errorbars are 
inversely proportional to the significance of the peak. As expected, the position of 
the low-lying peak increases as the average valence quark mass is increased. The 
dotted vertical lines passing through each of these peaks indicate the values for 
the ground state energies obtained from double cosh fits to the correlators [29] . 
Therefore, for the groundstates we find that the results of the MEM analysis are 

6 Meson spectral functions using domain wall fermions in quenched QCD have been studied 
in Ref. [28]. 
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Figure 7: Pseudoscalar spectral functions in QCD with 2 + 1 flavours of dynam- 
ical domain wall fermions, determined using the Maximum Entropy Method, for 
different values of the valence quark masses. The vertical dotted lines indicate 
groundstate masses obtained with conventional cosh fits [29|. The inset shows a 
blow-up of the second bump. 



in precise agreement with the results of conventional fitting techniques. 

At higher frequencies, 1 < auj < 2, second much broader bumps are visible. 
There appears to be no significant valence quark mass dependence of the position 
of these second peaks. We note that the heights of the bumps are of the same 
order as the structure due to the cusps observed in the free fermion calculation, for 
a domain wall height of unity0 From our earlier analysis, we find therefore that 
the position of these bumps is not inconsistent with their identification as lattice 
artefacts. However, these peaks may also contain excited state contributions, 
but due to the width of the peaks these are difficult to resolve. An unambiguous 
way to disentangle excited state resonances from the contribution due to lattice 
artefacts is to carry out an analysis at different lattice spacings. With current 
data, this option is not yet availablejf] 

7 Note, however, that we have not applied wave function renormalization, which affects the 
vertical scale. Moreover, the free calculation indicates a strong dependence on the domain 
wall height, such that a quantitative comparison would require knowledge of the renormalized 
domain wall height. 

8 Unphysical structure at higher frequencies has also been observed in quenched simulations 
with Wilson fermions. In Rcfs. 42 , 413] this was interpreted as bound states of Wilson doublers. 
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7 Summary 



We have analyzed meson spectral functions from lattice fermions with chiral 
symmetry. For free fermions, we have given a general prescription on how to 
construct lattice meson spectral functions from the euclidean fermion propaga- 
tor, extending the analysis of Ref. [I]. We have subsequently applied this to 
overlap fermions, domain wall fermions and overlap hypercube fermions. Lattice 
artefacts appear at (twice the) frequencies at which the lattice dispersion rela- 
tion Ek deviates from the continuum relation. The most pronounced effect is 
the appearance of cusps at frequencies determined by dE^/dh = 0. For most 
fermion formulations (Wilson, staggered, standard overlap, domain wall), these 
cusps appear at frequencies 1 < au < 2. In order to shift these artefacts to higher 
energies, it is necessary to use lattice fermions with improved scaling behaviour, 
such as hypercube fermions. In our spectral function analysis, we found that 
using the hypercube operator as a kernel in the overlap formalism indeed yields a 
formulation with good chiral and scaling behaviour, as could be anticipated from 
previous studies [TT] . 

From a comparison between overlap and domain wall spectral functions, we 
found that the latter have a remarkably strong dependence on the domain wall 
height, whereas the dependence on the corresponding parameter in the case of 
the overlap operator is much milder. For a domain wall height of unity and a 
small fermion mass, we found that the overlap and domain wall spectral functions 
are comparable. 

In the final section of the paper we have performed a Maximum Entropy 
analysis of pseudoscalar spectral functions in QCD with dynamical domain wall 
fermions, using data generated by the UKQCD and RBC collaborations. We 
found good agreement between the groundstate masses, determined by conven- 
tional cosh fits, and the position of the peak in the spectral functions. At larger 
frequencies, 1 < au < 2, a second peak can be seen. We have argued that 
this structure is not inconsistent with the lattice artefacts found in the analytical 
study, although the presence of excited states cannot be excluded. An unambigu- 
ous way to distinguish (physical) excited states from (unphysical) lattice artefacts 
discussed here, is to repeat the analysis at different lattice spacings. 
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A More on overlap hypercube fermions 

In this appendix we collect some expressions for overlap hypercube fermions, 
discussed in Section El 

The pole is determined by 

C* F \K) sin 2 k v = -m 2 b 2 (K) } b(K) < 0, (A.l) 

V 

where m was defined in Eq. (13.71) . In order to solve for the dispersion relation, 
we follow Ref. [3] and use the notation 

Cf F (K) = 5k, 
C i F (K) sin hi = ctjk + (3^ cos k^, 

m H F(-^) = ^ik + ^2k cos k 4} (A. 2) 



a i*> ^ = H "k " Ac = ^ a ikAk- (A.3) 



and 

3 

' 9 ^,9 , 

«k 

i=l i=l i=l 

Expressions for a^, /3; k , <5k, and Ki^k can easily be derived from these definitions 
combined with Eqs. (I5.3[ 15.41) . Explicit expressions are given in Eqs. (B.1-B.9) 
of Ref. [3]. 

Writing /c 4 = iE k yields again a quadratic equation for coshi^k, where the 
allowed solution is of the form 



with 



cosh^k = — , (A.4) 

^k 



x k = «k • Ac + "^o K 2k («ik - A 1 ) , 

l/k = K-/3k) 2 +(^-^) (Sl + oi) 

+rhl [2 (a k • Ac) ^2k («ik - At) + («ik - A 1 ) 2 ( 5 k ~ © ~ «L ( s l + a Q] > 
^k = ^ - Pi - m\t^.. (A.5) 

The root (xk — -^/l/kV-Zk < for all momenta inside the Brillouin zone and 
therefore not a valid solution, in contrast to the standard HF case. 
The rest mass is given by 



m^o («io - /-O + 5 J 51 + ml \(k w - fj) 



\ 2 _ .2 
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coshE = — ' (A-6) 

"o — m K 2O 

22 



where 



Kl0 = A (0) + 6A« + + 8A^, (A.7) 

K20 = 2A (1) + 12A {2) + 24A (3) + 16A (4) . (A.8) 



The explicit expression for the residue reads 

1 fi 1 w(zJ5k, k) 



(A.9) 



£ k /i 2 + mo/4 - cosh E k - a k ■ /3 k - ml^biiE^ k) sinh E k 

There exists no pole at fc 4 = n — iE k , provided b(K) = mwp{K) — fi > 0, or 

A 4 < K ik — ^2k cosh E k . (A. 10) 

Since «i,2k depend on in a nontrivial manner on the coefficients in the HF kernel, 
this constraint on [i has to be verified on a case by case situation. 
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